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Abstract 

The dynamics of tilted vortex lines in Josephson-coupled layered supercon- 
ductors is considered within the time-dependent Ginzburg-Landau theory. 
The frequency and angular dependences of the complex-valued vortex mobil- 
ity \x are studied. The components of the viscosity and inertial mass tensors 
are found to increase essentially for magnetic field orientations close to the lay- 
ers. For superconducting/normal metal multilayers the frequency (uS) range 
is shown to exist where the value depends logarithmically on lo. 

PACS numbers: 74.80.Dm, 74.20.De, 74.60.Ge 
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The inertial mass M and viscosity rj are known to be important quantities for char- 
acterizing vortex dynamic response in type-II superconductors. Without taking account of 
pinning effects, the vortex line bending and traction by the superflow the simple phenomeno- 
logical description of the vortex dynamics in isotropic superconductors (see, e.g. Ref. [Q) 
is based on the following equation for the flux line velocity V: MV + rjV = <pojtr x n/c, 
where j tr = jo^xp(iujt) is the transport current density, n is the unit vector which points 
in the magnetic field direction, 0o is the flux quantum. In the context of many applica- 
tions (including high-frequency phenomena such as microwave or infrared response) it is 
convenient to introduce also the complex-valued dynamic mobility fi = (iujM + rj)^ 1 which 
may be extended to include both pinning and flux-creep effects 0,0]. For T close to T c the 
quantity rj may be estimated within the Bardeen-Stephen model or calculated more exactly 
using the time-dependent Ginzburg-Landau (TDGL) theory H. The inertial effects may 
play an essential role not only in high-frequency classical dynamics but also in quantum 
tunneling phenomena at extremely low temperatures |J. There exist several contributions 
to the vortex inertial mass |5Hl2|, i.e., one from the electronic states within the normal core 
(the electronic contribution M c ), one connected with the electric field energy (M em ) and 
one arising from mechanical stress and strain effects. Note that for isotropic compounds the 
term M em is much less than the electronic contribution M c ~ h 2 Nf (Nf is the density of 
states at the Fermi level). 

Recently a great deal of attention has been devoted to the peculiarities of the vortex 
dynamics in extremely anisotropic high-T c superconductors which pertain to a larger class 
of materials including quasi-two dimensional (2D) organic superconductors, chalcogenides 
and artificial superconducting multilayers. A common feature of these systems is the weak 
interlayer Josephson coupling which results in a short effective coherence length £ c for the 
order parameter spatial variation along the c direction (perpendicular to the layers). In a 
broad temperature range the £ c value may be much smaller than the interlayer distance D. 
This fact results in the quasi-2D character of static and dynamic magnetic properties in 
these compounds. As a consequence, the description of the dissipative and inertial effects in 
the vortex motion within the anisotropic TDGL theory employing an effective mass tensor 
fl3HT6| appears to be adequate only over a limited temperature range (£ C (T) > D). In this 
paper we restrict ourselves to the case of weak magnetic fields when the so-called effective 
cores of neighbouring flux lines do not overlap (see Refs. [|17],[il|) and study the motion of an 
isolated vortex line which forms the angle 7 with the z axis (c direction) and lies, e.g., in the 
(xz) plane. It is the purpose of this Letter to report on a calculation of the vortex dynamic 
mobility for layered superconducting structures in the temperature range £ C (T) <C D. Let us 
consider a stack of thin superconducting (S) films of thickness d separated by insulating (I) 
or normal metal (N) layers of thickness D ^> d. The inertial effects and a viscous drag force 
acting on a moving flux line are determined by the complicated structure of the effective 
vortex core RI7| , |I8| . For tawy < £/D the effective core size is of the order of the coherence 
length £ in the plane (xy). In the angular domain Dtanj > £ the effective core consists of 
2D core regions connected by the Josephson-like vortex cores which appear for Dtan'j > Lj , 
where Lj = (Dhc 2 / (&ne\ 2 ab J c )) l l 2 is the Josephson length, J c is the interlayer Josephson 
critical current density, X a b is the London penetration depth for currents along layers. The 
characteristic dimension of the 2D core region in the plane (xy) is = min^Dtawy, Lj]. 
The dimensions of the Josephson-like core along y and z are Lj and D, respectively. 
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Before presenting the derivation of the vortex dynamic equation, we briefly discuss the 
main mechanisms responsible for the inertial and dissipative effects in the vortex motion. 
One may expect that the contributions arising from the regions of 2D pancake normal cores 
do not change essentially and may be calculated in a manner analogous to the treatment 
in Refs. p|-|T0||. The novel specific terms in M and rj come from the 2D and Josephson-like 



core regions. In particular, the simple qualitative arguments show that with an increase in 
7 the term M em should increase due to the capacitive effects. This term is determined by 
the spatial distribution of the electric field generated by a moving flux line. Far away from 
the normal core regions and for rather large distances between the neighbouring 2D pancake 
vortices (Dtawy > £) the averaged z component of the electric field E zn between layers n 
and ri + 1 may be found using the Josephson relation (we take here the gauge A z = 0): 
E zn = —%W{9 n — n+ i)/ (2eD), where 9 n (r — R(i)) is the phase of the complex order pa- 
rameter in the layer n, R = V is the vortex velocity. At large distances R n from the 2D 
vortex center (in the 2D core region) E zn decreases as R~ x to the distances R n ~ a^. Such an 
extremely slow decay of E zn results in the logarithmic divergence of the kinetic energy of the 
electric field, which is cut off at a large length scale ad- Thus the contribution proportional to 
the value In (ad/£) should appear in the expression for the inertial mass. This enhancement 
of the M em value may be significant if the dielectric constant of insulating layers is large. 
The analogous logarithmic term from the 2D core region contributes to the viscosity coeffi- 
cient due to the dissipation produced by the interlayer normal currents (see also Ref. [|l9f| ) . 
Surely, the influence of these interlayer currents on the scalar electrodynamic potential dis- 
tribution (and, hence, on the complex- valued dynamic mobility) is most significant for S/N 
multilayers and temperatures close to T c . For large angles tawy > Lj/D the formation of 
Josephson-like cores results in the additional contributions to the vortex viscosity and mass, 
which should be linear in the length Dtawy. These terms may play an essential role in the 
vortex dynamics as it follows from the results obtained in Refs. |||20| for the particular case 
7 = 7r/2. Obviously, the Josephson strings parallel to the x axis contribute to the vortex 
mobility only for a nonzero velocity component V y . As a consequence, for large tilting angles 
the vortex mass and viscosity coefficient should depend on the transport current orientation 
in the plane (xy). 

We now continue with the calculation of the viscosity coefficient and the vortex mass 



within the phenomenological Lawrence- Doniach model |2T| generalized for the description 
of the dynamic phenomena (for T close to T c ). Such a generalization may be written in 
analogy to the TDGL theory in 3D homogeneous superconductors. In the following treat- 
ment we will use the diffusion type TDGL equations which are strictly valid only for gapless 
superconductors. Nevertheless, we believe the final results to be qualitatively correct for a 
more general case. Note also that in the TDGL equations given below we neglect the terms 
connected with the particle-hole asymmetry and responsible for the vortex traction by the 
superflow [^,^3|]. If we choose the gauge A z = 0, the equation for the order parameter 
ip n (r) = f n exp(i9 n ) and the current density (averaged over the periodicity length D) have 
the form: 

/ 2ie \ 2 

= £ 2 ( v - — a„j ij n + ij n 
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- \^n\ 2 ^n + J^i^n+l + VVi-1 ~ ^n) (1) 

hc 2 f 2 2e 1 • 

jn|1 = 8^~ (W " ~ h~c A ^ ~ aab{S7Vn + c An) ' (2) 

Here (p n is the electrochemical potential in the n-th layer, <j ah is the normal state conductivity 
in the layer direction, r = 7ih/(8(T c — T)). The interlayer current density is given by the 
expression: 

(jz)n,n+l = Jcfnfn+lSin(8 n+hn ) - O C D^ (f > n+1 >n (3) 

where a c is the interlayer normal state conductivity, 9 n +i )n = n+ i —9 n , <p n +x,n — fn+i ~Vn- 
Employing the expressions for S/N multilayers we assume that all the relevant length 

scales in the plane (xy) are larger than the periodicity length D. As a consequence, our 
approach is strictly valid only in the limit £ 3> D [pl| . We consider the dynamics of a tilted 



vortex line in the presence of an applied transport current J n (£) flowing along the plane (xy). 
For simplicity we suggest that the total stack thickness is much smaller than the effective skin 
depth and the current density J n = J ir does not depend on the layer number n. Considering 
only straight tilted vortices, we also leave aside the problem of possible vortex line flexures 
immediately under the surface. According to Ref. |p4| , for the particular case Dtawy Lj 



these flexures occur in the boundary layer of the width L s ~ rran[A a &, a]D/Lj <^ X a b, where 
a is the intervortex spacing. The generalization of the technique discussed below, taking 
account of the vortex line bending is straightforward (see, e.g., Ref. |fl). 

To obtain the equation of motion we use the procedure similar to the one discussed in 
Refs. |3],|l9|,^] and search for the solution of Eq.(|l|) in the form: 



i ! n = G n + g n = G(r - R - r n ) +g(r-R- r„, t), (4) 

where r n = (nDtgj, 0), G(r — r„) is the order parameter distribution for a static vortex line 
and g is the small correction which is of the order of V = R. The arguments analogous 
to the ones used for isotropic type-II superconductors (4| show that in the limit \ a b 3> D 
the magnetic field generated by a moving vortex line may be neglected within the effective 
core region (which provides the main contribution to the inertial mass and viscosity). The 
condition X ab ^> D may be easily met in real layered structures, at least, for temperatures 
close to T c . As a consequence, in Eqs. we may take into account only the vector 

potential A n = — 47rA^c _1 j s generated by the transport supercurrent density j s averaged 
over the structure period. Substituting the expansion @ into Eq. ([!]) and neglecting the 
higher order terms in V one obtains the equation for the correction g n . We next multiply 
this equation by G pn = pVG n (p is an arbitrary translation vector in the plane (xy)) and 
integrate over the area S of the circle |r — r n — R| < Ri where R\ meets the condition 
R\ 3> Dtawy. After simple transformations the real part of the resulting equation reads: 



— P(J* 

c 



x z ) = ^-Re J G* pn (G vn - y^ n G n - g n )d 2 r, (5) 



where r] = 0.5a a bU(p H C 2C 2 , H c2 is the upper critical field for 7 = 0, 
u = he 2 1 (?>2X 2 ab o a b(T c — T)), and G vn = VVG„. The numerical factor u is determined by 
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the pair-breaking mechanism (see Refs. [p5| , p6[ for details). The term in the left-hand side 
(l.h.s.) of Eq. (ID corresponds to the Lorentz force acting on a 2D pancake vortex while the 
right-hand side (r.h.s.) contains the terms responsible for the viscosity and inertial mass. 
We do not consider the extrinsic forces due to pinning and the interaction with other vortex 
lines, which add to the l.h.s. of Eq. (|5]). Let us examine at first the contribution to the 
integral in the r.h.s. of Eq. @ which comes from the domain |r — R — r„| ^> £ (outside the 
normal cores). In this case one can put \ip n \ ~ 1. The continuity equation for the layered 
system reads: 

dPn j. . \jz)n,n+l (jz)n— l,n n , c \ 

-Qjr + dtvj n + = 0, (6) 

Pn = -^^{2^ n -^ n+1 -<p n -i), (7) 

where p n is the averaged charge density in the layer n and e is the high-frequency dielectric 
constant. We neglect here the difference between the electrochemical and scalar electrody- 
namic potentials, i.e., assume the Thomas- Fermi screening length to be less than all the 
relevant length scales. Using Eqs. ([[]), (§P,(0) one obtains: 

— $„ = £ 2 A6» n + j^(sin9 n+1>n - sin9 n ^i) (8) 

$n = \S+ J^rJ^J (Vn+l + Vn-l ~ tyn), (9) 

where s = cr c £ 2 / (ua a bD 2 ) and $ n = ip n + H9 n / (2e) is the gauge-invariant scalar potential. 
Note that we assumed u > 1 which is consistent with the microscopic theory results for gap- 
less superconductors p5| , ^6| . Let us consider the two-term expansion for the order parameter 
phase: 6 n = 6 nv + \n = d v (r — R — r n ) + Xm where 9 v (r — r n ) is the phase distribution for 
a static vortex line and Xn is the first-order in V correction. To solve the system ©,© 
we use here the linear approximation for the interlayer Josephson current density and re- 
place the terms sin6 n+ i^ n by 9 n+ i >n . The validity of such an approximation was discussed 
in detail in Ref. JTS|. The solution of the system (|H),(]9) gives us a possibility to evaluate 



the integral in the r.h.s. of Eq. (y) far away from the normal cores taking account of the 
expressions G pn ~ «(pV 9 nv )exp(i9 nv ) , g n ~ ix n exp(i9 nv ) . Note that the contribution which 
comes from the normal core region (|r — r n — R| ^ £) does not differ essentially from the one 
obtained previously for a vortex line in isotropic superconductors Finally one obtains 

the following equation of vortex line motion in the Fourier representation: 

/i^ 1 V;(^)x + ii y l V y {u)yQ = ^-3 tr (uj) x z , (10) 

where xo, yo, zo are the unit vectors of the coordinate system. We follow here the treatment 
in Refs. [|l],[7| and include both the supercurrent and normal current densities to the driving 
Lorentz force in the r.h.s. of Eq. (0). The complex- valued dynamic mobilities n XyV appear 
to depend strongly on the parameter s and tilting angle 7. The s parameter may be written 
in the form: s = l^/D 2 , where the length l c = C,\Jcr c /ui7ab is the electric field penetration 
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depth along z. For S/I multilayers one has l c ~ £ c and, as a consequence, s < 1 in the 
temperature range considered in this paper. For this limit one obtains: 



fx- 1 = [tcu(M c + M 2D )+r lc + r l2D ]D- 1 (11) 
a y 1 = I 1 ® 1 + (ivMj + r]j)tanj (12) 

M 2D = r-rprMl + a d /0 ; r/ 2D = (13) 

where ?7 C = tjqociD, M c = r]oTa 2 D and ai, a 2 are the constants of the order unity. The 
terms w c and iuM c are connected with the dissipation and inertial effects in normal core 
domains. The logarithmic terms r] 2 o, M 2 £> in Eqs. (13) may be considered as the con- 



tributions to the viscosity and inertial mass of a single 2D pancake vortex. Let us com- 
pare the M 2 £> value to the inertial mass M c of the normal core for high-T c materials (for 
the angular domain tawy 3> £/D). Considering Bi — 2 : 2 : 1 : 2 as an example we take 
D ~ 15A, e~10, v~b {T ~ T c ) ~ 10 _4 ficm, £(T = 0) = £ - 20 - 40A, T c ~ 80 - 90K, 
Lj/D ~ 300 - 1000 and obtain 

— — ~ — — - u — /n — ~ 0.1 In — 

M c r n l a 2 huD i o b £ £ 

Thus, for tawy > Lj/D one has M 2 £, ~ M c . This rough estimate suggests that the con- 
tribution M 2 £, to the vortex mass of electromagnetic origin may be significant in layered 
high-T c superconductors. In the angular domain Dtawy ^> Lj the vortex dynamic mobility 
depends essentially on the orientation of the current with respect to the in-plane magnetic 
field component. Such a dependence results from the contribution to /i" 1 , proportional to 
the length Dtan r y of the Josephson string connecting the neighbouring 2D pancake vortices. 
The terms rjj and Mj correspond to the viscosity and inertial mass per unit length of the 



Josephson vortex. Note that the expressions (14) are in good agreement with the ones 
obtained previously in Refs. PJ2H| for vortices parallel to the layers. 

The normal currents in nonsuperconducting layers lead to the penetration of the field E 
generated by a 2D vortex (moving in the plane z = nD) at a finite length l c along z. For 
S/N multilayers the l c value may be larger than the structure period D. This fact results 
in an essential decrease of the potential ip n in the 2D core regions (£ < |r — R — r n | < ad). 
As a consequence, the order parameter phase in these regions satisfies the diffusion type 
equation (see Eq. (§)), where Dg = £ 2 /t is the diffusion constant. The characteristic time 
scale of the phase distortion propagation through the 2D core region is to ~ a \lDe- If the 
frequency uj of the applied ac field is higher than the value u ~ t$ , then the essential time 
dispersion comes into play and the logarithmic divergence of the mobility coefficient is cut 
off at a length scale L u determined by the uj value. The upper cutoff is of the order of 



a diffusion length L w = ^Dq/uj. The validity of these qualitative arguments is proved by 
direct calculations. As an example we consider here only the limit l c ^> D (s ^> 1) which is 
relevant to S/N multilayers, at least, for temperatures close to T c . We also restrict ourselves 
to the frequency range ujt <^ 1 (L w ^> £). The complex- valued dynamic mobility has the 
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form //~* = (l + iu)T/u)r] Xj y, where the quantities M X)V = tj] XjV /u and 7] XtV may be considered 
as the components of the effective mass and viscosity tensors, respectively. Evaluating the 
integral in Eq. (|5|) to the logarithmic accuracy one obtains: 

Vx ~r] (F + Q/q); Vy ~ T] (F + qQ) (15) 



1 + (Dtan-f/Lj) 2 - F = ft + ln(l + L m /£) (16) 
/ min[L WJ a d ^/s]\ 

Q = ln { 1+ £ + Lm )+&> 

where |/3i j2 | Z 1, L m = mm[a d ,Lj. The logarithmic in u terms contribute to the dy- 
namic mobility as well as to the effective mass and viscosity in the frequency range 
s _1 ^ 2 a^ 2 < cur < 1. For ^fsDtan^ < £ the logarithmic terms in Eqs. flTBD,(|r7D are small 
and the effective mass and viscosity are determined by the dynamic processes in the normal 
core regions. In the opposite limit (^Dtaw-f > £) these terms are large and the values 
Pi, 02 may be neglected in a wide frequency range. The latter conclusion is not valid only 
for Dtawy 3> Lj and L u Lj, when the logarithmic term in Eq. (|17|) vanishes while the 
value 02 ~ —iL^/Lj provides an essential imaginary contribution to the coefficient rj y . As 
a consequence, the M y and rj y values introduced above become complex-valued and, hence, 
in this case it is not adequate to define these quantities as the effective mass and viscosity 
(otherwise these definitions may be used). For higher frequences lot > 1 the vortex mobility 
is completely determined by the order parameter dynamics in the normal core region, which 
was analysed, e.g., in Ref. 0. 

In summary, we have obtained the equation of tilted vortex motion, which may be appli- 
cable to various layered systems, including high-T c copper oxides. The angular dependences 
of the viscosity coefficient and the inertial mass have been investigated. For S/N mutilayers 
the complicated effective core structure is shown to result in the specific frequency depen- 
dence of the complex- valued dynamic mobility even in the frequency range out < 1. The 
approach developed above provides a starting point for the study of high-frequency linear 
response in the mixed state of layered superconductors. For this purpose the vortex dynamic 
mobility should be certainly extended to include both pinning and flux-creep effects 
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